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Correcting second-order contamination in low-resolution spectra 
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An empirical method for correcting low-resolution astronomical spectra for second-order contamination is presented. The 
method was developed for correcting spectra obtained with grism#4 of the ALFOSC spectrograph at the Nordic Optical 
Telescope and the performance is demonstrated on spectra of two nearby bright Type la supernovae. 
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1 Introduction 

The availability of large format CCDs with high quantum 
efficiency over a wide wavelength range, makes possible 
most modern low-resolution spectrographs to have the ca- 
pability to obtain spectra that cover the whole optical range 
(3200-10000A) in a single exposure. However, it follows 
from the theory of diffraction gratings (see, e.g.,Schroeder 
2000) that different diffraction orders overlap, i.e. a photon 
with wavelength A m in the m-th order will be diffracted at 
the same direction as a photon with wavelength A m+1 from 
the m-|- 1-st order and thus both will be recorded at the same 
pixel on the detector. For diffraction gratings the relation 
between A m and A m+1 is simple, A m = (to + l)\ m+1 /m 
(but see Gutierrez-Moreno et al. (1994) for a case where this 
does not hold true due the a specific spectrograph design). 
Many spectrographs employ grisms instead of gratings, in 
which case the wavelength overlap relation is not that sim- 
ple and is generally a non-linear function: 



A' 



Most low-resolution astronomical spectra are obtained in 
the first diffraction order (with typical dispersion ~ 2 — 
5 Apixel -1 ) and thus beyond ~ 6000 A will be contami- 
nated by the blue light of the 2nd-order, in many cases quite 
significantly. The order contamination greatly compromises 
the spectrophotometric accuracy but also affects the mea- 
sured equivalent widths and profiles of the lines, as well as 
may introduce spurious lines when emission line objects are 
studied (Gutierrez-Moreno et al. 1994). The traditional way 
to overcome this problem has been to use blue light block- 
ing filters or two different spectrometer settings. In either 
case one either loses the blue part of the spectrum or doubles 
the exposure time. Double-beam spectrographs are another 
solution, however, these are not as common instruments as 
the single-beam spectrographs. Besides, the dichroics used 
to split the beam may have non-uniform response at cer- 
tain wavelengths, making the accurate flux calibration chal- 
lenging. It would be therefore beneficial in many cases if 



a methodology for correcting single spectra for the higher- 
order contamination is available. However, little has been 
done in that direction so far. In previous works Bowers et 
al. (1997), Peroux et al. (2001) and Norman at al. (2002) 
only briefly outline the procedures they use to correct the 
order contamination in their spectra, and Gutierrez-Moreno 
et al. (1994), Leonard et al. (2002) and Steinmetz et al. 
(2006) discuss other aspects of the order contamination. I 
started working on the development of a method for 2nd- 
order correction in early 2004, largely motivated by the de- 
sire to optimize the Target-of-Opportunity observations of 
Type la Supernovae (SNe la) that began at the Nordic Opti- 
cal Telescope (NOT) in the fall of 2003 using the ALFOSC 
spectrograph. Because grism #4 (3200-9100 A) suffers of 
severe 2nd-order contamination beyond ~ 5800 A, obser- 
vations with the red grism #5 (5000-9800 A; free of order 
contamination) were also obtained, but this doubled the ob- 
serving time. Eventually, in the end of 2004 Szokoly et al. 
(2004) published a paper where they also presented in de- 
tail a method for correcting for the 2nd-order contamina- 
tion; but my work was done independently, albeit arriving 
at quite a similar method. 

2 The method 

The observed spectrum in units of detected counts, A^A 1 ), 
is N{\ 1 ) = A^A 1 ) + A II (A I ), where A^A 1 ) is the true 
lst-order spectrum that would be detected if there was no 
contamination and A 11 (A 1 ) is the 2nd-order contamination, 
which is the true 2nd-order spectrum A 11 (A 11 ) mapped on 
the lst-order wavelengths A 1 with Eq.Q] The efficiency of 
a grisrrQ at a given order should depend only on the pa- 
rameters of the grim (Palmer & Loewen 2006fl In this 



1 Throughout the paper I will refer to grisms, but the discussions are 
valid for diffraction grating as well. 

2 The efficiency of reflection gratings depend sensitively on the po- 
larization of the incidence light. However, the transmission gratings, to 
which class grisms belong, are practically free from the polarization effects 
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Fig. 1 He arc-lamp spectra a) with and b) without the 
blocking filter GG475. c) The combined spectrum obtained 
through the U and B filters. The lines used to obtain the 
wavelength solution for the second order are labeled in the 
c) panel. Some of these lines are clearly seen in the spec- 
trum without the blocking filter. 



case it is easy to show that for any object observed, the 
2nd-order photons with wavelength A 11 that are recorded at 
wavelength A 1 in the first order, 7Y II (A I ), can be expressed 
as a function of the detected 1 st-order photons N (X s -): 

N Il (X l ) = C , (A I )/V I (A I ). (2) 

The function C(A I ) is the ratio of the efficiencies^ of 
the 2nd- and lst-order C(X I ) = E n (X I )/E I (X I ), where 
E II (X I ) is the efficiency of the 2nd-order again mapped 
on the lst-order wavelengths. Thus the application of the 
method involves to steps: 1) wavelength transformation of 
the observed spectrum according to Eq.Q] i.e. N(X l ) — > 
7Y(/(n->i)(,\ 1 )) and 2) multiplying 7Y(/f tI_>I )(A 1 )) by 
C(A I ) and subtracting it from A^A 1 ). What needs to be de- 
termined is the wavelength relation Eq.Q]and C(A r ). 

The work on the method was first initiated by the ob- 
servations of SN2004S (Krisciunas et al. 2007) that were 
obtained at NOT in February 2004 using ALFOSC plus 
grism #4. Because of the ~ —31° declination, the obser- 
vations at La Palma were performed at very high airmass, 



(Palmer & Loewen 2006). Besides, in most cases the observed objects will 
not be strongly polarized. 

3 efficiency here means the fraction of the incident photons in the cor- 
responding order that are recorded on the detector. 
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Fig. 2 The overlap relation between the 1 st and 2nd order. 
The lcr scatter around the fit is 0.2 A. Inset: the non-linear 
part of the fit. 



often X > 2, with slit along the parallactic angle to min- 
imize the slit losses due to the differential atmospheric re- 
fraction (Filippenko 1982). At such a high airmass, the at- 
mospheric refraction is large and the atmosphere effectively 
acts as a cross-disperser, similar to cross-dispersers in the 
Echelle spectrographs. This allowed on one occasion the 
2nd-order to be spatially separated from the first order (see 
Fig. [51 top), and the spectra at the two orders could be in- 
dividually extracted. My first attempt was to use this obser- 
vation to derive a method for correcting for the contami- 
nation. With too few blue photons detected (because of the 
high airmass) this was unsuccessful, but stimulated further 
work using observations of bright, blue stars with and with- 
out order-blocking filters in order to isolate the 2nd-order 
light. N u (X l ) is the difference between the photon flux de- 
tected without and with the filter. Such observations were 
only obtained in November 2005. Using a slightly adjusted 
theoretical wavelength overlap function for grism #4 pro- 
vided by Per Rasmussen (Copenhagen Observatory, private 
communication) the method was successfully tested. How- 
ever, more observations were obtained in May 2006, in- 
cluding dedicated observations to determine experimentally 
the wavelength overlap relation Eq.Q] The results presented 
here are based on this later observing run. 

To determine the overlap relation Eq.Q] observations of 
Ne and He arc-lamps were obtained. The lst-order wave- 
length solution was obtained from exposure of the He and 
Ne lamps simultaneously. Four 3 min long exposures of the 
He lamp without any filter and through the U, B, and the fil- 
ter GG475 that block the light shortward of ~ 4500 A were 
also obtained. The U and B filters served to block the lst- 
order He lines, so that the faint 2nd-order lines cold be mea- 
sured to derive the wavelength solution for the 2nd-order 
(Fig. [I])- The other two exposures were only used to check 
if the filters introduced shifts to the line positions. Shifts of 
1-2 pixels were found (and corrected) for the U and B fil- 
ters, and none for GG475. Figure[2] shows the positions of 
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the 2nd-order He lines as detected in the lst-order vs. their 
real wavelengths. These data determine the order overlap 
relation and as can be seen from the inset in Fig.|2]it is non- 
linear. Note also that because of the higher dispersion, more 
lines could be detected in the second order than are usually 
seen in the first order. 

Figure[3h shows the observations of the spectrophoto- 
metric standard Feige 66 with and without the GG475 fil- 
ter. Figure^ shows the difference between the two spec- 
tra (= ^"(A 1 )^ and the wavelength transformed observed 
spectrum A^A 1 ) -» ^(/("^(A 1 )). Note that the Hydro- 
gen B aimer lines from the 2nd-order are clearly seen in the 
difference spectrum and their position are well matched by 
the lines in the transformed spectrum. Figure[3j; shows the 
ratio N u (X I )/N(f ( - u ^ I ' ) (X I ) which is the correction func- 
tion C(X l ). Because of the low sensitivity in the blue, the 
part between 5800A and 6000A is rather noisy and I fit a 
parametrized function (the thick line), which is the final cor- 
rection function C(X l ). 

3 Performance 

To demonstrate the performance of the algorithm I use ob- 
servations of two Type la supernovae (SNe) 2005cf (Gar- 
avini et al. 2007) and 2005hk (Stanishev et al. 2007). Both 
SNe were observed at NOT with grisms #4 and #5 as a part 
of large observing campaigns involving many other tele- 
scopes. If there were no order contamination in grism #4, 
the shape of the flux calibrated spectra with the two grisms 
would be the same. I therefore correct the grism #4 spectra 
for the 2nd-order contamination and compare them with the 
ones obtained with grism #5 in order to evaluate the per- 
formance of the algorithm. Note that the observing condi- 
tions were probably not perfectly photometric and so there 
is difference in the absolute flux level of the grism #4 and 
#5 spectra of the order of ~ 5 — 15%. Such small varia- 
tions should not afect the relative flux calibration and are 
not problem for the analysis presented here. 

Figured presents the results for SN 2005cf. In Fig.@h 
are shown the observed spectrum with grism #4, the es- 
timated correction and the corrected spectrum. Figure^ 
shows the uncorrected, flux calibrated spectra with the two 
grisms. Note that for this example neither the SN nor the 
spectrum of the spectrophotometric standard BD +33 2642 
used for the flux calibration were corrected. Because 
BD +33 2642 is actually bluer than the SN, effectively flux 
is subtracted from and not added to the SN spectrum. After 
correcting for the 2nd-order contamination, the spectra with 
the two grisms line up nicely (Fig. [4}:). Figure^ clearly 
demonstrates not only the problem with the accurate relative 



4 The GG475 filter transmission decreases nearly linearly from 97% to 
92% between 5800 A and 9000 A. The flux observed through the blocking 
filter was corrected accordingly. This correction was not in the original 
method and was included later after a note in Szokoly et al. (2004) that 
the blue light blocking filter may also affect the red part of the spectrum. 
This is the only influence that the paper of Szokoly et al. had on the work 
presented here. 
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Fig. 3 a) The spectra of the hot sdO star Feige 66 with 
and without the GG475 blocking filter, b) The 2nd-order 
light and the wavelength transformed lst-order light, c) The 
ratio between the 2nd and 1 st-order light. The parametrized 
fit which is used for the 2nd-order contamination correction 
is also shown. 



flux calibration, but also the negative effect on the measured 
line equivalent widths: in the uncorrected grism #4 spectrum 
the strong Ca II infrared triplet at ~ 8000A is much weaker 
than in the grism #5 spectrum. 

With respect to SN 2005hk, I first correct the spec- 
trophotometric standards to derive a contamination-free re- 
sponse function, which is then used to flux calibrate the SN 
spectra. On 2005 October 31 SN 2005hk was observed at 
high airmass and similarly to the case of the SN 2004S ob- 
servations described in the previous section, the second or- 
der can be seen spatially well separated from the first or- 
der (Fig. [5] top). The three ticks in Fig. [5] (top) roughly cor- 
respond to wavelengths of ^6700, 7900, and 9000 A from 
left to right. Cuts along the spatial axis were taken at these 
three positions in order to measure the FWHM of the spec- 
tral traces and the separation between them. The separa- 
tions are 12, 10 and 7 pixels from blue to red wavelengths 
(1 pixel=0.19 arcsec). The FWHM of the lst-order trace is 
4.6 pixies at all cuts, but the second order trace FWHM de- 
creases from 6 to 5 pixels from blue to red. The large sep- 
aration between the orders allows to extract the lst-order 
light with almost no 2nd-order contamination and to com- 
pare how the correction works not only with the grism #5 
spectrum but also with this 'clean' lst-order extraction. I 
used an extraction aperture of 9 pixels width and thus be- 
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Fig. 4 a) The observed, the calculated correction and the 
corrected spectra of SN 2005cf observed at NOT on 2005 
June 2.b) The flux calibrated uncorrected spectrum with 
grism #4 compared with the spectrum with grism #5. The 
effect of the 2nd-order contamination in grism #4 is very 
strong, c) Comparison of the two spectra after the one with 
grism #4 has been corrected for the 2nd-order contamina- 
tion. The two spectra match almost perfectly. 

low ~ 8000 A the extracted spectrum should be to a large 
extent contamination-free, while above ~ 8000 A the con- 
taminations will be somehow larger because the separation 
between the traces decreases to the red. Nonetheless, I find 
that the so extracted grism #4 spectrum is in a good agree- 
ment with the grism #5 one, indicating that nearly 'clean' 
extraction have been achieved. The results are presented in 
Fig. [5] and turning our attention to the 2nd-order contamina- 
tion correction, I again find that the contamination has been 
very well corrected. 

It should be noted that the experience from the observa- 
tions shows that separation between the orders as large as 10 
pixels is very rare. It is normally of the order of few pixels 
only and the two orders are blended togetheJI. For the algo- 
rithm to work correctly, the spectrum that is to be corrected 
has to contain the entire light from both, the 1 st and the 2nd 
order. Given the possibility for a significant separation be- 
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Fig. 5 Upper panel: The two-dimensional spectral image 
of SN 2005hk observed at NOT on 2005 October 31 with 
grism #4. The 1st and 2nd orders are clearly seen spatially 
separated (see the text for details), a) The flux calibrated un- 
corrected spectrum with grism #4 compared with the spec- 
trum with grism #5. b) Comparison of the two spectra after 
the one with grism #4 has been corrected for the 2nd-order 
contamination, c) Comparison between the grism #5 spec- 
trum with the uncorrected grism #4 spectrum. In this case a 
narrow extraction aperture was used for grism #4 spectrum, 
so that only the 1st order light was extracted. 



tween the orders, to make sure that this is indeed the case 
one should carefully select the aperture position and width 
so that the light in both orders is fully extracted. 

4 Conclusions 

I present an empirical method for correcting low-resolution 
astronomical spectra for 2nd-order contamination. The 
method was applied to grism #4 of the ALFOSC instrument 
at the Nordic Optical Telescope, but can be used at any low- 
resolution grating/grism spectrograph. The performance of 
the method was tested on observations of two bright nearby 
Type la supernovae, SNe 2005cf and 2005hk, and in both 
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cases the result was excellent, greatly improving the accu- 
racy of the final flux calibrated spectra. Thus, in many sit- 
uations the method would allow to reduce almost twice the 
observing time of programs that need to observe the whole 
optical range. Various Target-of-Opportunity programs on 
Supernovae, Gamma Ray Bursts and other transient object 
can greatly benefit from it. I expect the method to work very 
well for low-polarized objects with broad spectral features 
like SNe and GRBs. Its performance on objects with narrow 
spectral features and/or high polarization should be further 
investigated. 
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